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Abstract 

It is demonstrated that quasi-exactly solvable models of quantum mechanics admit 
an interesting duality transformation which changes the form of their potentials and 
inverts the sign of all the exactly calculable energy levels. This transformation helps one 
to reveal some new features of quasi-exactly solvable models and associated orthogonal 
polynomials. 
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1 Introduction 



In this paper we demonstrate that quasi-exactly solvable models of quantum mechanics (i.e. 
models admitting exact solutions only for certain limited parts of the spectrum)^ normally 
appear in the form of dublets. The models forming a given dublet look differently but their 
exactly calculable energy levels exactly coincide up to a sign transformation. We call such 
models dual to each other. Note that the duality property does not hold for other (exactly 
noncalculable) levels. 

In this paper, we consider simplest quasi-exactly solvable models introduced in ref. and 
represented by even polynomials of degree six (section 2). We construct the dual pairs of these 
models (section 3) and show that the duality property enables one to estimate the intervals 
of the spectra occupied by exactly calculable energy levels (section 4). We also demonstrate 
that the duality property can be used for finding classical counterparts of the phenomenon 
of quasi-exact solvability (section 5) and for establishing new remarkable properties of non- 
standard orthogonal polynomials associated with quasi-exactly solvable models (section 6). 
A brief discussion of the general case will be given in the concluding section 7. 

2 The model 

The Schrodinger equation for the simplest sequence of quasi-exactly solvable models reads 

j-Zi^^ + x2(ax2 + 6)2 - ha{2M + 3)x^^ ^'(x) = ^'f(x), (2.1) 

where b is real, a > is a positive real number and M = 0, 1, . . . is an arbitrary non-negative 
integer. 

The explicit solutions of this equation can be represented in the form 



E = h 
and 

M 



M 

1 



(2M + l)6 + 2a^^, 



(2.2) 



i=l 



1 / hx^ ax^ 



(2.3) 



where the numbers i = 1, . . . , M (playing the role of wavefunction zeros) satisfy the system 
of numerical equations 



— = 6ei + <f, i = l,...,M (2.4) 

k=i 



and the following additional condition 

M 



Y.(^ = 0. (2.5) 



i=l 



For details see original papers [Q, |2| ^ and reviews ^. 
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The proof of formulas (|2.2| ) - (|2.5| ) is trivial. For details see the proof of analogous formulas 
in ref. @. The number of exactly calculable energy levels and wavefunctions is, obviously, 
given by the number of solutions of the system ( |2.4| ) - ( |2.5D . At first sight, this system is 
overdetermined, because it contains M + 1 equations for only M unknowns ^j, i = 1, . . . , M. 
However, it can be shown that, due to the symmetry of this system under reflections .^j — > — ^j, 
it still has non-trivial solutions for any M, and their number is equal to [M/2] + 1, where 
[ ] is a standard notation for an integer part. Note that the parity of the obtained energy 
levels coincides with the parity of M. It can also be shown that the numbers satisfying 
the system ( |2.4] ) - ( |2.5| ) are either real or purely imaginary. 

3 Duality and self-duality 

It is easily seen that formulas (|2.1| ) - ( p. 51) for the "plus" model (with b = \b\) can be reduced 
to the same formulas for "minus" model (with b = —\b\) after the changes 

X ix, i^i, E -E. (3.1) 

Hereafter we shall call the models (^]^) , differing from each other by the sign of the parameter 
6, dual to each other and the transformation ( ^.Ip — the duality transformation. We also 
shall use the notations 

V{x) = x^{ax^ + bf - ha{2M + 3)x^ (3.2) 

for the potential of the original model ( |2.lD , and 

V{x) = x^{ax^ - bf - ha{2M + 3)x^ (3.3) 

for its dual. If the parameter b is zero then the potentials V{x) and V{x) exactly coincide 
and read 

V{x) = V{x) = a^x^ - ha{2M + 3)x'^. (3.4) 
It is natural to call the model ( |3.4[ ) self-dual. 

4 Spectral properties 

From formula (^]^) it follows that, if E is an exactly calculable energy level in model (|3.2| ), 
then E = —E is also a certain exactly calculable energy level in the dual model (p.3|). It 
is interesting to establish a relationship between the quantum numbers (ordinal numbers) 
of these levels. Assume that the energy level E describes the Kth excitation, E = Ex- 
Then, according to the oscillation theorem, the wavefunction (|2.3| ), ^(x) = should 
have K real zeros (nodes) and, therefore, K numbers should lie on the real x-axis. In 
this case, the M — K remaining numbers should lie on the imaginary x-axis. But from 
formula (^]^) it follows that, after the duality transformation, all the imaginary numbers 
describing the imaginary wavefunction zeros in the model ( |3.2| ) become real, = i^j, and 
take the meaning of real wavefunction zeros (nodes) in the dual model (|3.3| ) . But this means 
that ^{x) = ^ m-k{x) the corresponding energy level describes the (M — K)ih. excitation, 
E = Em-k- This leads us to the simple relations 

Ek = -Em-k, "^Kix) ~ ^M-K{ix). (4.1) 
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In particular, the minimal exactly calculable energy level in a given quasi-exactly solvable 
model ( |3.2D coincides with the maximal exactly calculable energy level in the corresponding 
dual model ( |3.3D taken with sign minus! For the self-dual model ( p. 41 ) the situation becomes 
even more transparent. In this case, because of the coincidence of the sets {E} and {E}, we 
obtain 

Ek = -Em-k, ^k{x) ~ '^A4-K{ix). (4.2) 

The abovementioned facts enable one to detemine the intervals of the spectra of models 
(|3.2| ), (|3.3| ) and occupied by the exactly calculable levels. We denote these intervals by 
[Emin,Emax] and [Emin,Emax], respectively. The duality property means that 

E-max — EYnirii EjYiax — EYnin- (^•^) 

On the other hand, we have the obvious inequalities 

Emin > min V{x), Emin > min V{x). (4.4) 
Comparing ( [4.3[ ) with ([4.4|), we obtain 

Emax < -min V{x), E^ax < -min V{x). (4.5) 

This means that 

[Emin, Emax] C [min V{x), -min V{x)], [Emin, Emax] C [min V{x), -min V{x)], (4.6) 
where 

{0, b > ^/ac, 

2(27a)-i (-2h + V62 + 3ca) (-3ca + 6^ + foV^^ + 3ca) , h < ^^'^^ 

and 

{0, —h > ^/ac^ 

2(27a)-i(+26 + \/6^T3^) (-3ca + 62_5V6^T3^), -6 < ^^'^^ 

Here we introduced a new quantity 

c = (2M + 3)7l. (4.9) 



For the self-dual model ( |3.4| ), formulas ( |4.6| ) - ( f4.8[ ) become simpler. Taking in them 6 = 0, 
we obtain 



[Emin, Emax] — [Emin, Emax] C 



(4.10) 



The examples of the dual and self-dual potentials with the intervals occupied by exactly 
calculable energy levels are depicted in figures 1 and 2. 
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Figure 1: The form of the potential V{x) = x^iai? + 5)^ — acx^ (left picture) and its dual 
V{x) = x^{ax^ — b)'^ — acx^ (right picture) for a = 1, 6 = 3.3 and c = 3.5. The gray 
domains correspond to the parts of spectra occupied by the exactly calculable levels. Here 

[Emin, Emax] = [0,12.37] and [Emin, Emax] = [—12.37,0]. 

5 The classical limit 

It is naturally to ask what is the classical analogue of the phenomenon of quasi-exact solvabil- 
ity of models (p.l|)? As far as we know, this question has never been raised in the literature. 
To answer it, we should understand what is the difference between the solutions of the classical 
equation of motion corresponding to the cases with E S [Emin, Emax] and E [Emin, Emax]- 
For this it is very helpful to use the duality property. 

First note that for small h the spacing between the neighbouring exactly calculable energy 
levels is of order of ^. If ^ — > then the lowest exactly calculable levels in models ( |3.2[) and 
( |3.3| ) approach the absolute minima of potentials V{x) and V{x). For this reason, the classical 
limit of the relation ( [4.6| ) reads 

[Emin,Emax] = [min Vo{x), -min Vb(x)], [Emin,Em.ax] = [min 1^0(2;), -min ^0(2;)], (5.1) 

where by Vo(x) and Vo{x) we denoted the classical limits of the potentials V{x) and V{x). 

For determining the classical potentials remember that the original ones contain two 
arbitrary parameters H and M. Therefore, the result of classical limit will depend of how 
fast M tends to infinity if h tends to zero. This forces one to consider separaely two cases: 
1) ?i ^ and c ^ and 2) h ^ and c 7^ 0. These limits we shall call respectively the soft 
and hard classical limits. 



The soft classical limit. In the soft classical limit the potentials (3^) and (|3.3| ) take 
the form 

Vo{x) = x^{ax^ + bf, Vo{x) = x^{ax^ -bf. (5.2) 

In this case we have 

min Vo{x) = 0, min Vo{x) = 0, (5-3) 

and therefore, the parts of the spectrum occupied by the exactly calculable levels are collaps- 
ing into a single (ground state) energy level E = E = 0. For this reason, it is naturally to 
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Figure 2: The form of the self-dual potential V{x) = V{x) = a^x^ — acx"^ for a = 1 and 
c = 3.5. The gray domain corresponds to the part of the spectrum occupied by exactly 
calculable levels. Here [Emm, Emax] = [—2.52, +2.52]. 



expect that the classical analogue of the phenomenon of quasi-exact solvability should lie in 
a comparative simplicity of the ground state solution {E = E = 0) for the model (|3.2|), ( |3.3| ) 
with respect to all other solutions {E > 0). 

In order to make sure that this is really so, consider, for example, the expression for the 
subbarier classical action S{x) in the model (^.2|). Remember, that function S{x) is a natural 
classical counterpart of the quantum wavefunction ^(x) and is related to it by the formula 
^{x) ~ e~^^^^^^. Solving the corresponding Hamilton - Jacobi equation we obtain 

S{x) = J dx^x^iax^ + b)^ - E, E>0. (5.4) 

In some sense this is an explicit expression which, after evaluating the integral, can be reduced 
to the combination of tabulated elliptic integrals. Let us however consider the particular case 
with E = when some analogue of quasi-exact solvability is expected. In this case the 
integral ( ^.41) can be evaluated explicitly and the result takes a very elegant and simple form 

5(x) = ^ + ^, E = 0. (5.5) 

We see that this case is distinguished from all other cases by the fact that the action function 
is in this case regular, without any branch - points and cuts! Comparing the special solution 



( |5.5| ) with the general one (5.4), we can say with confidence that the former has all rights to 
be called "exact solution". 

The hard classical limit. In the hard classical limit the potentials of the models ( |3.2| ) 
and ( ^.3| ) take the form 

Voix) = x'^iax'^ + bf - acx^, Voix) = x'^{ax'^ - bf - acx^ . (5.6) 

The number of exactly calculable energy levels tends to infinity and the spacing between them 
tends to zero. The sets of former quantum exactly calculable energy levels transform into 



the segments of a continuous spectrum given by the formula (5T). The subbarier classical 
action now reads 

S{x) = / dx^J x'^{ax'^ + 6)^ — acx^ — E. (5.7) 
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It is not difficult to see that, in the complex x-plane, the function S{x) has six square- 
root-type branch point singularities corresponding to the roots of the sixth-order algebraic 
equation x'^{ax'^ + 6)^ — cax"^ = E. These singularities are connected pairwise by cuts. A 
simple analysis shows that if £J £ [Emim Emax], then all the six singularities lie on the real and 
imaginary x-axes. However, ii E ^ [Emin, Emax], then there are singularities lying outside 
the axes. The location of these singularities and the corresponding cuts are depicted in the 
figure 3. 
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Figure 3: The complex rr-plane for the classical action S{x) (formula ( |5.7D ). Here, as before, 
a = 1, b = 3.3, c = 3.5 and [Emin, Emax] = [0, 12.37] is an interval occupied by the former 
exactly calculable energy levels. In the left picture E = 1 £ [Emin, Emax], while in the right 
picture E = 20 ^ [Emin, Emax]- We see, that the disctribution of singularities and cuts for 
these two cases is different. 

So, we see that there is some qualitative difference between the solutions of classical 
equations corresponding to the former exactly calculable and exactly non-calculable energy 
levels. At the same time, there is no difference between the complexity of these solutions. 
This is quite obvious, because if M — > oo the algebraic equations determining the exactly 
calculable states in quantum quasi-exactly solvable models become infinitely complicated and 
the difference between them and any other states dissapears. 



6 The duaUty and orthogonal polynomials 



The wavefunctions corresponding to exactly calculable energy levels in models (|3.2| ) and (|3 
can be represented in the form 



b ± rr 







1 




exp 


~h 






1 




exp 





ix^ bx"^ 



(6.1) 



where p = 0, 1 is the parity of the solution, and Pm{t), m = 0, 1, . . . , [M/2] are certain 
polynomials of degree [M/2]. From the condition of the orthogonality of wavefunctions 

f + OO 



'^m{x)'^n{x)dx = 0, 



^m{x)^n{x)dx = 0, m ^ n 



(6.2) 
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it follows the orthogonality of polynomials Pm{i) ctnd Pm{t): 

/ Pm{t)Pnit)L0{t)dt = 0, / P^{t)Pn{t)u{t)dt = 0, m^Tl (6.3) 

JO Jo 
with the weight functions 

uj{t) = tP-he-4-b\ oj{t) = tP-h~^+^\ (6.4) 



Note now that the relation (4J) for wavefunctions imphes the analogous relation for the 
corresponding polynomials 

Pm{t) ~ pM-mi-t) (6.5) 

whose comparison with (|6.3D gives 



oo rO 

Pm{t)Pn{tMt)dt = 0, / Pm{t)Pn(t)uj{t)dt = 0, 171 ^ n. (6.6) 
J-oo 

So, we see that polynomials Pm{x), m = 0, 1, . . . , [M/2] form an orthogonal system with the 
weight function lv{x) on two different intervals [— oo,0] and [0, oo]. 

Note that the existence of two intervals on which the polynomials Pn{x), n = 0, . . . , [M/2] 
form an orthogonal system is closely related to the facts that all these polynomials are of the 
same degree [M/2] and their number [M/2] + 1 is finite. Indeed, in this case, the total number 
of coefficients of these polynomials (excluding the leading ones which can be considered as 
normalization factors) is [M/2] ([M/2] + 1). Requiring the orthogonality of polynomials on a 
certain interval, we obtain [M/2] ([M/2] + l)/2 equations for their coefficients. In this case 
the number of equations is only one half of the total number of unknowns which prevents 
one from determining these polynomials uniquely. However, if we require the orthogonality 
of these polynomials on two different intervals, then the number of equations will exactly 
coincide with the number of unknowns, and the polynomials can be determined uniquely. 
Note that this property is typical for all non-standard orthogonal polynomials associated 
with quasi-exactly solvable models from ref. [^. 

7 Conclusion 

The duality property of the simplest quasi-exactly solvable model ( ^j]) and its applications 
discussed in sections 3-6 can easily be extended to all one-dimensional quasi-exactly solvable 
models listed in ref. ||3|. In order to demonstrate this, remember that these models are 
associated with the class of equations 

Loat + d)t{t -a)-^ + {At"^ + Bt + C)—- M{M + ^ - 1) ^ P{t) = EP{t) (7.1) 

(or their degenerate forms) having M -|- 1 polynomial solutions for any M. The transition of 
equation ( |7. 1| ) to the Schrodinger form can be performed after introducing new variables x 
and new functions ^'(x) by the formulas 

^= [ (72) 
J -uj{t + a)t{t - a) 
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and 



w N r.. . fdx\^ i f At'^ + Bt + C ) , 

in which the variable t belongs to the interval [0, a] (for u; > 0) or to the interval [a, oo] (for 
a; < 0). The set of parameters to, a, d, A, B, C determines the form of the resulting potential. 
From (7.1) it immediately follows that the transformation t — > — t, E —E reduces the quasi- 



exactly solvable model characterized by the set of parameters {u;,a,d,A,B,C} to another 
model characterized by the set {uj, d, a, A, —B, C}. This transformation is nothing else than 
the generalization of the duality transformation for model (^.l|). 

It would be an interesting problem to investigate Lie algebraic origin of the duality prop- 
erty of quasi-exactly solvable models and find a multi-dimensional generalization of this phe- 
nomenon. We intend to consider these questions in fortcoming publications. 
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